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Abstract 

By use of single instanton approximation the twist-3 light-cone distribution amplitudes of pion 
as well as the p-wave mesons, i.e., /o(1370), Kq(1430) and ao(1450) are investigated within the 
framework of QCD moment sum rules with inclusion of instanton effects based on the valence quark 
model. Results show that there is much change on light-cone distribution amplitudes by the isospin- 
and chirality-dependent instanton contribution compared with the instanton-free ones. We find the 
intanton involved twist-3 LCDAs are non-positive-definite within some range of momentum fraction 
and there are rapid changes at two ends of momentum fraction. To guarantee the convergence of 
moments by the method in this work a low instanton density should be adopted, for instance tIc = 
^fm~^ is welcomed. Possible ingredients which might have impact on the results are briefly discussed. 
These light-cone distribution amplitudes may be helpful to analyze exclusive heavy flavored processes. 

PACS numbers: 
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I. INTRODUCTION 



The hadronic light-cone distribution amphtude(LCDA or hght-cone wave functions) which 
parameterize the non-perturbative effects play an important role in understanding exclusive 
hard processes in QCDjliJai. The LCDA is also one of the key ingredients in QCD factoriza- 



tion approach 



J] for describing exclusive hadronic B decays which are helpful to ascertain the 



quark content of light scalars 



-7j] and other observables of phenomenological interest. There- 
fore if there is more information on LCDA there will be more complete study on heavy-flavored 
exclusive processes consequently many important parameters in standard model. To extract 
LCDA the nonperturbative method should be employed due to its nonperturbative nature. The 



QCD sum rules|8|] have been proved to be very successful one in obtaining useful information 
in hadronic region {q]. While it was found later after its born in order to produce reasonable 
results of lowest-lying states of pseudoscalar channel the instanton effect 10HlJ| should be taken 
into account in QCD sum rules 151]^. Then the role of instanton in QCD sum rules is exten- 



sively investigated 
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. Recently it was noticed that instanton may be helpful in lifting the 
mass degeneration in light scalar mesons above lGeVtl9|. The calculation of LCDAs within 

n 

the framework of QCD sum rules was introduced in[20| by studying the correlation function of 
currents with derivatives in this way the desired Gegenbauer coefficients can be expressed in 
terms of the moments derived from QCD sum rules. Using this method Chernyak and Zhit- 
nitsky(CZ) were able to describe many experimental data available at that timej2l|. Marr ying 
with lattice simulation the first two Gegenbauer coefficients of pion and kaon were studied 
under this method. 

In analyzing B decays to p-wave mesons 



22| 



23 



241] (more complete review on B decays, see 25|) 



such as /o(1370), -^"0(1430), ao(1450) the LCDAs are important input. Having noticed the 
significant irnplication of instanton in producing the realistic mass spectral of /o, Kq and 



above IGeV 



in this paper we will make some effort to investigate the instanton effects on 



the twist-three LCDAs of these mesons by the CZ method. The work in this paper can be 



^ For thorough review on instanton in QCD one can refer to[l6|. 
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regarded as partly pursuing of Ref. [l9| and the quark content are assigned as follows 



fo = -^{uu + dd), 
Kq = ds, 

ao = -^{uu-dd), (1) 

The LCDAs of pion, as a widely used quantity in exclusive processes, naturally attach attentions 
based on various methods from instanton model 261]. To this point we would like to stress that 



instanton is crucial in reproducing reasonable results of pion and its partners in 0~ channel 
by QCD sum rules[15|. Perturbatively the intuition is that due to the chirality suppress there 
should be larger size of ss component in u than rj. If it were the case u would be close to 
but the fact is that u nearly degenerates with p while there is considerable mass gap between 
u and 0. The empirical spectrum is that t] is nearly pure SU(3) octet while much smaller ss 
component in u. Then there arise the puzzle why there is much larger splitting between pion 
and rj than that of p and u which is difficult to explain perturbatively. The puzzle can be 
solved if the instanton is considered. Due to the definite chirality of fermion zero mode there 
is direct instanton contribution to pseudoscalar channel but not in the vector and tensor ones 
and additionally this instanton-induced correction is isospin-dependent thus pion and rj are 
lighter than p and u and there is larger mass gap between pion and t] than p and u. Hence 
pion provides a natural laboratory to study instanton and it is expected there should be some 
effects of instanton on its twist-3 LCDAs. Therefore it is meaningful to analyze the twist-3 pion 
LCDAs by CZ method with inclusion of instanton. In fact the isospin-dependence of instanton 
effects is also important to reproduce a realistic mass gap between an and /o in 0+ channel 



above IGeV by QCD sum rules assuming small size ss mixing into fo 19|. The quark content 
of pion is as usual 

tt'^ = —= [uu — dd) . 

It is hoped these twist-3 LCDAs may shed some light on the exclusive hadronic B decays. We 
would like to emphasize that in QCD language a real hadron should be described by a set of 
Fock states which each one has the same quantum number as the hadron. For example 

\Ko) = i^glds) + ^g,\d-sg) + iJ^,^^\d-sqq) + •-, (2) 
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It is no doubt there are also twist-3 LCDAs introduced by higher Fock states but we will not 
consider them here. In other words we deal with leading Fock states or the valence quarks of 
the hadrons. Additionally there were some efforts in obtaining pion LCDAs from holographic 



QCD methods 



271 ] and nonlocal condensates based on QCD sum rules 



as well as some work 



on the shape of pion LCDAs 29|]. These work gives some insights of pion LCDAs to us. 

The outline of this paper is as follows: in Sec. II we present the sum rules with inclusion of 
instanton contributions to obtain the LCDAs of /o, Kq, oq and pion. In Sec. Ill the numerical 
results and discussions will be presented, finally we summarize our conclusions in Sec. IV. An 
appendix is given to show the vanishing of instanton contribution to tensor moment sum rules 
within the method. 



II. BASIC FORMULAS 

A. moments and light-cone distribution amplitudes 

Firstly we define the decay constants of scalar meson S and the pseudoscalar meson P 

{^\Mi\S{p)) = msfs, 
{0\q2^1,ql\P{p)) " ^'''"'^ 



mi + m2 

where mg, mp mi and ma are the mass of scalar meson, pseudoscalar meson, qi and ga, 
respectively. The twist-3 LCDAs (psi'^) ^"^^ for the scalar meson S with quark content 

giga are defined asp] 

{0\q2{z2)qi{zi)\S{p)) = msfs fdue'^^^^-^^+^^-^^^Un), (3) 



{0\q2{z2WMzi)\S{p)) = -msfs{p,z,-PuZ,) jyue'^^^-'^^"-^-^^^^^, (4) 
the two twist-3 two-particle LCDAs (f)^p{u) and (f)p{u) of pseudoscalar meson are defined as [3^ 
{0\q2{z2)^l,ql{zl)\P{p)) = ['due'^^'^-^^+^^-'^UW, (5) 



mi + ma Jq 



.2 pi 



{0\q2{z2Ku76qi{zi)\P{p)) = -\^^^{p,z,-PuZ,) I rfwe^("^-^^+^^-^^V^(«), (6) 

3 mi ma Jo 
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where u always refers to the momentum fraction carried by one quark and u = 1 — u is another 
quark momentum fraction; z = Z2 — z\. Noticing the gauge- invariant Wilson path-ordered 
integral 

[z2,zi] = Pexp[ig / da^Ai^ia)]. 

J Zl 

has been suppressed. The normalization of these four twist-3 LCDAs are 

1 .1 
du(j)g{u) = / du(f)g{u) = 1, 
Jo 

du(tfp{u) = I du(t)p{u) = 1, (7) 
Jo 

To proceed firstly we remind the reader that there is no pure zero-mode contribution to tensor 
moment sum rules(see the appendix) since instanton effect is chirality-dependent thus we do 
not consider (psiu) and (ppiu) here. For the instanton-free tensor moment sum rules to calculate 
twist-3 LCDAs of p-wave mesons above IGeV one can refer to Ref. 3l|. For simplicity in this 
paper we study the two 01 (m) and ^^(m). In fact we can only concentrate on the scalar moment 
sum rules because the pseudoscalar one can be deduced from the scalar one by some appropriate 
substitutions, therefore in following we deal with the scalar sum rules only. 
Generally the twist-three LCDA <f)s{u) has the following form 

00 

0|(n,/i) = 1 + ^a„(/i)Cy2(2n - 1), (8) 

n=l 

1 /2 

where CJ {x) is Gegenbauer polynomials of order 1/2, the lowest ones are [32 1 

Cl/\x) = ^x'-h, Cl/'(x) = i(35x^-30x2 + 3), (9) 



and the orthogonality relation is 



C'J'{x)CK'{x)dx = -^6nm. (10) 

2n + 1 



From Eq.Q one can easily derive 



mmtz-1^rq2ms{p)) = msfs{p-zr{c), (n) 
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where 



iC) 



rfM(2M-l)>|(M,/i). 



(12) 



From the orthogonal relation Eq.l lTOl) the Gegenbauer moments a„ can be expressed in terms 
of (^"), for our purpose 



a, = ^-{3{e) - 1), a, = ^(35(a - 30(e^) + 3). 



(13) 



The next step is to calculate the so-called moments appearing in Eq. flTT]) . to this end we consider 
the following two-point correlation function with derivatives 



(14) 



Inoiz^q) = I J rfW^^(0|TO„(x)Ot(0)|0) 

= iz-qriTiq') 



with 



On{x) = qi{x){iz ■ D rq^ix), 0^(0) = g2(0)gi(0), 



The above correlation function can be expressed in terms of the operator product expansion. 
Up to leading order of and dimension-six we get^ 



-g^ In 



q^ 3 + n .a 



+ 



1 _ 1 _ 

m2{gsqi(rGqi) - —^mi{gsq2(rGq2) 



+ 



2q 

4:71 a. 



2g4 



27 

48 a, ,^ 
9"q 



'^ri2 + 3n-4)((gigi)2+(g2g2)') 



— -i{QiQi){Q2q2) 



(15) 



with n is even thus only the scalar even moments exist. This is the theoretical side of the 
correlation function from the quark-gluon dynamics. On the other hand Eq. fll4p can also be 



^ Noticing the operator product expansion is different from Ref. 31| on the mass-dependent condensates terms, 
but there is httle impact on the results since these terms are greatly suppressed by the quark mass. We find 
the operator product expansion in Eq. (|15p presents a well extremum behavior. 
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derived phenomenologically based on the dispersion relation 

1 /"°° lmP''(s) 

lJq^) = - ds ^^ + subtr. const., (16) 

TT Jo s - g2 

The imaginary part Im/^ol-^) is obtained by inserting a complete quantum sets ^ \'n){n) into 
Eq.f ll4p which reads 

Im4l(g^) = vrm|/|(r)5(g^ - m|) + V^(9^ - ^o), (17) 

By equating the theoretical and phenomenological sides of Ino{z, q) we get the sum rules 

Substituting Eq.f llSp and Eq. lfTTll into Eq.( IT8l) . taking Borel transformation and subtracting 
the continuum contributions we arrive the desired scalar moment sum rules 

— + (gigi) + (^mi H ^"^2 j (g2g2) 



2M2 2M2 
487r 

^]^:^Wi^i/W2g2). (19) 

Up to now all our analysis still confines in the conventional QCD sum rules, in the coming 
subsection the instanton will take part in the game. 

B. Inclusion of the instanton effects in moment sum rules 

In this subsection our calculation are in four- dimension Euclidean space unless explicitly 
point out. Instanton is the nontrivial solution of classical field equation in four-dimension 



Euclidean gauge-field theories which is first discovered by Belavin et al 10|. Subsequently 't 
Hooft 3| derived the instanton with topological charge Q = 1 in four-dimension Euclidean 
space 

A» = -Va,. . (20) 

^ g [x — xo) + p 
G%{x) = --Va,u- -2, (21) 



where p is instanton size, rja^y is the t'Hooft 77 symbol, Xq is an any point in four-dimension 
Euchdean space called instanton center. In this instanton background field there is quark 
zero-mode which represent the tunneling effects, for our purpose we write the quark zero-mode 
propagator explicitly in regular gauge^ 

.2 1 1 



P 



[{x-xoy + p'f [{y-xoy + p'] 



3/2 



7/.7i^^(l -75) ®[r+rj. 



where m* is the effective mass and 



(22) 



T,, 



{r,±^), 



(23) 



with the useful relations 



r r 



(24) 



The instanton contribution to the scalar moment sum rules is obtained by substituting the 
zero-mode propagator Eq.f l22|) into correlation function Eq.(IT] 



d^xe'^''{0\TOnix)O\0)\0) 



' [ d^xe'^^ I dpnip)p^ I d^xo^^ 



i2'l3 



with 



[x - XqY + p2 



0„(x) = qio{x){iz ■ '^rrq.oix), O^O) = g2o(0) 



3/2 



[x - Xo)^ + 



9l0, 



3/2 • 



(25) 



where the integration over collective coordinates of instanton is explicit, the anti- instanton 
contribution as well as traces over 7 and SU(2) matrix are completed implicitly. The instanton 
density has the simple form proposed in ISj 



n(p) = ncS{p - Pc), 



(26) 



For the propagator in singular gauge on can refer to [3; 



But now one should notice that the covariant derivative in Eq. fl2^ contains the instanton 
background field A^/^ to guarantee gauge invariance 



where is the SU(2) generator t"" = a"/2 with the normalization condition 

1 



(27) 



6' 



ab 



(28) 



Firstly we take a close look at the covariant derivative 

'3 1 - ^/) 



IZ 



(29) 



We can sort the terms in expansion of Eg. (12 9 P into two kinds. The first one only contains the 
differential operator which acts on zero- mode propagator from left and right, the second one 
is the remaining parts which include all the instanton field involved terms. In following we 
will elucidate the complete contribution from these two kinds at n = 2 and n = 4 by explicit 
calculation. 

Contribution to the correlation function of the first kind is easy to calculate such that we 
can derive a general formula for n as follows^ 



Ul''\Q^) = J rfW^^(0|TO„(x)Ot(0)|0) 



TT^ mtm* 



l<"'2 



[xl + p2)3 



X A{2iY 



[x^ + p 



2^n+3 ' 



(30) 



where A is constant 

^ = [i + (-i: 

n-l 



3 5 
2 ' 2 



2n + l 



I Y^^k^ ^y+k 3 5 2(n-A;) + l 3 5 2k + I 



k=l 



2 2 



2 2 



■* Without confusion we take pc ^ p for simplicity. 
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It is worthy to mention that in deriving Eq. fl30p from Eq. fl2^ the hght-cone constraint = 



is crucial otherwise there wil 
of the following formulae 32 



be big masses. Now the remaining work is trivial with the help 



34| 



^zdz 



z'^KJz] 



Tiu) V 2 ; 



v—m \^Ji 



2v-i 



KJz) = KJz) 



(31) 



where = —q^ and Ky{z) is the MacDonald function. When the smoke clears, we get the 
desired results for the first kind contribution 



n 



d^xe'^^<d\TOn{x)0\^)\^) 
p2 2"+i[l + (-1)"](1 +n) 



X {z-QrQ^Kl{Qp), 



n + l 



(32) 



We find the non-vanishing contribution of the second kind for n = 2 is 



]-[second^g2^ 



ciW^^5io(0, X] Xo){-2igz ■ A^^^fS2o{x, 0; 



(33) 



Substituting the zero-mode propagator Eq.( p2|) and instanton field Eq.( l20|) at regular gauge 
into Eq.f p3|) . for definite we write its explicit form as follows 

QAir^mlm^ J J J [Xq + p^Y'^ 

' 1 , , , _ 

It^lu- (1 - 75)r+r^ 



n 



second //O 2 



[{x - xqY + p"^]^'^ """2 



2 {X - Xo)5 2 [X - Xq)^ , 

^ ^cr VaaS , \ n . n ^ Zp Tjiypy \ 9 , 9 

g (x — xoj + p g [X — xq) + p 
1 1 1 , , , _ 

X ^^^^^ (l_^5)T-+r^ 



ncpt 



2 
1 



GAn^mlm*^ J J (x^ + plf [(a; _ ^0)2 + ^2] ■ 

X Vaas{x — XQ)sZprii,p^(^X — XQ).y 

X tr[7^7^^(l - 75)7a7/3^(l - 75)]tr [r+r~rVV+r^~] , 



(34) 
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After a lengthy calculation we obtain the contribution from instanton field involved part for 
n = 2 in expansion^ 

= 6^^" ■ Q^Q'^'^Qp)- (35) 

where the isospin dependence have been included. One can find that this part is comparable 
to the first kind contribution for n = 2. 

For n = 4 the situation is more complicated when we expand Eq.( l29l) order by order since 
the differential operator 9^ can also act on instanton field and give non-vanishing contribution. 
In considering this effects on instanton field we find the following terms for the second kind 
contribution to correlator 

nsecond^g2) _ ^4 ^ rf4^e^Q- 1 (Q, xo)4 ■ - ^)] '^20 (x, 0; Xq) | ( " 2z^7z ■ A^^)) ' 



-2igz ■ A^^^y 



S'2o(a;,0;a;o), (36) 



where for brevity the SU(2) generator index is suppressed. Some effort later we obtain 

nr^'HQ') = #^(^ ■ QrQ'KKQp). (37) 

2Umjm2 

Notice that we have included the anti-instanton effects both in Eq. fl5^ and Eq. fl57|) . Combining 
Eq.( l32p . Eq.f l35p and Eq. fl37p we get the complete zero- mode contribution to correlation function 

= nf^\Q') + -P^iz . QfQ'KliQp), (38) 

for n = 2 and 

urm = nr (Q^) + ^^^^(^ ■ QYQ'Kim. m 

for ?T, = 4, from the two equations above one can find the contribution of two kinds is comparable 
thus we can not omit one of them simply. 

After dealing with the quark zero-mode in the instanton background field, now we calculate 
non-zero mode contribution to the correlation function. For n = 2 this part is 

n^^"^(Q2) J rf4^e^'3^^™(0,x)(-2i^2 ■ A(^))2^f'^(x,0), (40) 



For simplicity we replace pc by p. 
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where S'f^™(0,x) and S'2^™(x,0) are the non-zero-mode propagators. In instanton background 
field the complete form of quark propagator consists of zero-mode and non-zero-mode parts 

Sj{x,y) = Sr{x,y) + Sr'^{x,y), (41) 

In previous paragraphs we have completed all contribution to correlation function from zero- 
mode in single instanton approximation. To obtain the non-zero-mode contribution we need 
S'j^™(x, y), while S'j^™(x, y) is a quite complicated object. It was shown|49l-l5l| that it is reliable 
to take massless free propagator approximation for Sf^'^{x,y) if the zero-mode contribution to 
Green function is maximal. Fortunately in present calculation this requirement can be met 
since there is direct instaton contribution to the correlation function we considered thus it is 
convenient to take massless free propagator approximation 

for non-zero-mode. At this point it is necessary to stress that in vector channel this approxi- 



mation is no longer valid otherwise the vector current is not always conserved 52| . In this case 



the propagator of non-zero-mode is very involved one can refer to Ref. 16|] for details. 
Combining Eq. lHOl) and Eq.( H2l) we arrive 



1 /■ ,4 [Z-{X- Xo)] 



2 



According to translation invariance {x — Xq) — )■ m we can factorize out the the integral with 
respect to instanton center xq thus there are two separate integral given by nonzero mode and 
instanton field. Since 2^ = it is easy to see the instanton integral is vanishing therefore there 
is no contribution of nonzero mode to correlation function. In other word nonzero mode and 
instanton field do not "entangle" each other. The physical meaning is that in this case the 
instanton field does not transfer momentum from one quark to another. Similar analysis also 
hold for n = 4. 

To be consistent with Eq.( IT9|) we should reformulate the total instanton induced contribu- 
tions in term of dispersion relation. For this purpose noticing the properties of MacDonald 



function under analytical continuation are 



35| 
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in above expression (z) is the Hankel function of the first kind 

Hl'\z) = Mz)+tY,{z). (45) 

where Ju{z) and Yiy{z) are the Bessel functions and Neumann functions, respectively. The last 
step is improving Eq.f l32p so that it is consistent with Eq.(IT9l). To this end noticing the cut 
structure of the Hankel functions Eq. (H5l) one can find 



lmK^{—ip^/s) = — Jl{p^/s)Yl{p^/s) + slugular term, (46) 

As usual in terms of the dispersion relation we get the final results improved by the Borel 
transformation and continuum contribution subtracted. To be definite we collect the whole 
results as follows 

,2 



mym)ew[-^] = nr^(.o,M^) 

^ ^ 3m* 



^/ rfs.Ji(pv^)yi(pv^)exp[-— ], (47) 

1 2 

for {C,s) and 

-l/l(aexp[-^] = nr^(.o,M^) 

+ ^-'^idkl '''M)y.(pVS)e.vl-j^], (48) 

for {^'^) where the isospin dependence has been considered and H^^^ is the right hand side 
of Eq.( IT9|) . It is easy to derive the pseudoscalar moment sum rules from the scalar ones if we 
consider the effects of ^75 carefully in calculations 

(mi + M^-' 



■^J\.sM,-fi)U,V-s)^.A-ipl (49) 



for {Q) and 



2 



"i^pfpi^p) r rrip-. jjps,oPE, ,^ 
(m, + m.)^ ^"P[~M^] ^^-^ 



ps,OPE 

,2 /-so 



4:5m 
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1^ / rfs. Ji(pv^)n(pv^) exp [ - — ] , (50) 



for {C,p) where 11^*''^^^ represents the following instanton-free sum rules 

3 1 r° r s , n-3,a 



nr 



OPE. 



^TT^ n + 1 Jq 
'n + 1 



24 ' TT 

n + 1 



I2q2) 



+ 



2M2 



m2{gsqicrGqi) 



;mi{gsq2crGq2) 



27 mA 
487r a, 

+ ——{qiqi}{q2q2}- 



2M2 



9 M2 

Obviously the instanton contribution in pseudoscalar channel is opposite to scalar one which 
reflects the chirality-dependence of instanton effects. 

Now all the formulae needed have been fixed. The parameters which will be adopted in our 



numerical analysis are as follows 



fly, 



37| 



a, = 0.517, (— G^) = 0.012 ±0.006GeV^ 

TT 

{uu) = {dd) = -(0.225 ± 0.15)^GeV^ {ss) = (0.8 ± 0.2){uu), 

rriu = 0.004GeV, = 0.006GeV, m, = 0.12GeV, 
sUaGu) = {gsdaGd) = 0.8GeV^(un), (ggSaGs) = 0.8{gsU(TGu) . 



(51) 



We read the masses and decay constants of /o, Kq and from 19 1 



ffo 


= 375MeV, 




= 370MeV, 


fao 


= 370MeV. 



(52) 



The mass and decay constant pion are 



37| 



= 140MeV, = 130MeV. 



(53) 



All the parameters in Eq. (l5T]) as well as Eq. fl52|) and Eq. fl53|) are taken at /i = IGeV. 

The remaining important parameters are the instanton related ones. For the effective masses 



14 



of quarks and instanton size the following values are work well 19|, |33| 



86MeV, 



P 



= 114±28MeV, 

-fm= — GeV~^ 
3 0.6 



The last parameter is the instanton density which still needs improvement. The original value 
Uc = |fni~^ is used widely 
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38 



40|. 



39| . while the lattice calculation suggested ric ~ lfm~'^ 
The work of Cristoforetti et a/j4l| based on the interacting instanton liquid model shown even 
a larger one was needed, i.e. ric = 3fm~^, to reproduce the nucleon mass and the low-energy 
constants in chiral perturbation theory. For this reason we will investigate the sensitivity of 
the moments consequently the LCDAs for different instanton density. 



III. RESULTS AND DISCUSSIONS 

Firstly we present the selection rule of the threshold and Borel window. The continuum(and 
exited states) as well as dimension-six condensates contribution should be controllable, as 
we demand that in the instanton-free sum rules the continuum contribution 



usual [8|, l31 



48 



the part in the dispersive integral from sq to oo should be less than 30% the total perturbative 
dispersion integration which sets an upper limit to us, the dimension-six condensates be no 
more than 15% which sets an lower limit to us. If there is extremum within the Borel window 
selected, we take it as our calculated value otherwise the mid-value within the window will 
be adopted. Then we turn on instanton contribution under same threshold and Borel window 
since in this way there will be well comparison for the two cases. Of course one can analyze 
the moments by separate threshold and Borel window when instanton effect turn on, however 
in general the sum rule is sensitive to threshold and Borel window that the instanton effect 
may be smeared by the change inducing by the new threshold and Borel window. In following 
discussions all the Borel windows satisfy the selection rule unless explicitly state. 

Along with the steps we find for pion the moment (^p) at threshold Sq = 4.0 ± 0.2GeV^ 
and Borel window e [1.35, 1.65] GeV^ as well as {0 at so = ± 0.2GeV2 and e 
[1.2, 1.5]GeV^ are stable from sum rules Eq.f l49l) and Eq.( l50|) at ric = respectively. The mid- 
value are (^p) = 0.34 and (^p) = 0.21 which are shown in Fig.[Tl When turning on the instanton 

15 



effects we find {C,p) = 0.52 and {C,p) = 0.36 at ric = ^fm~^. If increasing instanton density, for 
instance at ric = lfm~^ we find {^p) = 0.71 and {^p) = 0.50. Obviously the moments increase 
with instanton density increasing. Hence it is expected that at some larger Uc the second 



moment will be more than one, for instance if we take ric = 2fm~^ as proposed in 4l|] we find 
(^p) = 1.10 and (^p) = 0.80, it is astonishing that the second moment is more than 1! This 
result is unnatural since it is obvious from Eq.f lTS]) if assuming positive-definite LCDAs we have 

(^™) = f du{2u - l)"'(j){u, fi) < (C) = [ du{2u-lY(f){u,n) <1, (54) 
Jo Jo 

where 

n,m = 2, 4, , m > n. 

But the LCDAs themselves are not measurable, to observe their exact role we should convolute 
them with the hard scattering amplitudes T in exclusive processes. In fact the twist-2 LCDAs 



of /o(980) given in|5| from CZ method, twist-3 ones of pion 30| as well as the work in Ref. 42N46| 
show non-positive-definite behavior within some range of momentum fraction. In considering 
this it seems that our results may indicate non-positive-definite LCDAs. Indeed the instanton- 
involved LCDAs of pion show this property as presented in Fig.|3] although the moments are 
still well convergent. In considering the exact form of LCDAs and instanton density are not 
well known nowadays so that the impact of high instanton density to LCDAs by the CZ method 
seems to need further study. Our results show that it seems that the sum rules Eq. (H9l) and 
Eq.f l50p do not allow too large instanton density in order to get convergent moments. The 
obtained LCADs of pion for different instanton density is shown in the left panel of Fig.|3l 

Then we turn to analyze /o(1370) with the quark content assigned in Eq.([T]). The sum rules 
of this member is very similar to the pion except some different condensates terms induced 
by chirality. We can see the instanton contribution of /o(1370) is same as pion since there 
are simultaneous changes in its isospin and chirality relative to pion. At ric = we get the 
threshold and Borel window for (^J,) and {^j^) are sq = 4.7 ± 0.2GeV^, e [1.3, 1.7]GeV^ 
and So = 4.8 ± 0.2GeV^, e [1.8, 2.2]GeV^, respectively. Under the threshold and window 
there are well extremum behavior both for and (0^^) which are shown in Fig.[2j The 
extremum within the range of threshold are very stable we obtain (0^) = 0.35, (0^^) = 0.24. 
When the instanton effects turn on we find the mid-value are (0^^) = 0.55 and (0^) = 0.34 for 
ric = Ifiii""^ well as {0^) = 0.76 and (0^) = 0.44 for ric = Ifm"^. Similar to the case of 
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FIG. 1: Moments (^p)(left panel) and (^p) (right panel) of pion from instanton-free sum rules. Solid 
lines correspond to the central value of threshold while the dashed lines represent the threshold in- 
creasing by O.lGeV^ relative to central value(the same for Fig. [21 Fig. U] and Fig.[5|). 





1.95 2.00 2.05 

M'(GeV') 



FIG. 2: Moments (^^)(left panel) and (^f)(right panel) of /o(1370) at ^ = IGeV from instanton-free 
sum rules. 

71 at ric = 2fm~^ we find (.^|^) = 1.17 which more than 1. The moments also increase as the 
instanton density take a larger value which is understandable since the instanton contributions 
to pion and /o(1370) are equivalent due to the combined effects of isospin and chirality. Thus we 
conclude the instanton contribution is positive to the moment sum rules of pion and /o(1370). 
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FIG. 3: Twist-3 light-cone distribution amplitudes of pion(left panel) and /o(1370)(right panel) as 
function of momentum fraction u at = IGeV for different instanton density: solid line ric = 0, 
dashed ric = ^fm^'^ and dashed-dot = lfm~^, dash-dot-dot nc = 2fm~^(the same for Fig. [6]). 

The LCD As of /o(1370) are plotted in Fig.El 

It is easy to understand the similar impact of instanton on the LCDAs of pion and /o(980) 
since the dominant parts, i.e., perturbative dispersive integral and the instanton contribution 
of the moment sum rules for pion and /o(1370) are equivalent except the chirality-dependent 
condensates. In fact this similarity also reflects in the LCDAs which can be observed clearly 
from Fig.O One more important thing is that the LCDAs is positive-definite when there is 
no instanton effects while when the instanton involved there is strong impact on the profile 
of LCDAs. Due to the chirality-dependent parts at the two ends of momentum fraction the 
LCDAs of /o(1370) change more rapidly than pion. 

The sum rules of Kq[14:30) and ao(1450) are nearly the same since they share same isospin 
and chirality in addition to the difference introduced by flavor symmetry breaking. The adopted 
threshold and Borel window of (^|.,) and {^j^,) at Uc = are sq = 4.7 ± 0.2GeV^, e 
[1.3,1.7]GeV2 and sq = 5.5 ± 0.2GeV2, e [L35, 1.75]GeV2, respectively. There is little 
change of the extremum corresponding to different threshold both for the two moments, we 
get = 0.34 and (C.^*) = 0.23. When the instanton effects involved the moments change 

lot even at low density Uc = |fm~^. Both moments decrease compared with the case Uc = 0, 
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FIG. 4: Moments (^^)(left panel) and (^^) (right panel) of i^Q(1430) from instanton-free sum rules. 
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FIG. 5: Moments (^^)(left panel) and (i^g)(right panel) of ao(1450) from instanton-free sum rules at 
fi = IGeV. 



we obtain {^x*) = 0.17, {^k*) = 0.13. If the instanton density increases further, for instance 
Uc = lfm~^ we find there is flipping of the second and fourth moments, (^^.) = 0.01 and 
i^K^) ~ 0-04 which is unsatisfactory since it breaks the convergence. However at higher density 
rir = 2fm-^ the convergence can recover while it develops negative value - {^k*) = —0.32 and 
(^4,,) = -0.16. The LCDAs of K^{U30) for different instanton density are plotted in Fig.El 
It is clear that the profile of LCDAs of Kq{1A30) with ric ^ just reverse to pion and /o(1370) 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



0.0 0.1 0.2 0.3 0.4 0.6 0.6 0.7 0.8 0.9 1.0 „ 
u 

FIG. 6: Twist-3 light-cone distribution amplitudes of -fCg (1430) (left panel) and ao( 1450) (right panel) 
as function of momentum fraction u at /i = IGeV for different instanton density. 

which well indicates the conspiracy of chirality and isospin dependence of instanton effects. 

The case of ao(1450) runs in parallel to ^^0(1430). The threshold and working window 
of {^ao) ^^'^ (^ao) determined from with instanton-free sum rules are Sq = 4.9 ± 0.2GeV^, 
M2 g [1.5, 1.9]GeV2 and sq = 5.6 ± 0.2GeV2, A'P e [1.8, 2.2]GeV2, respectively. We obtain 
the extremum {C,^^) = 0.39 and (C^^^^) = 0.28 within the threshold range. After ric = |fni~^ 
turning on both moments are still well convergent and decrease to lower values: {^^o) ~ 0.20, 
(Cao) ~ 0.17. Similar to the case oiC,K*, at nc = lfm~^ we find there is also flipping of the second 
and fourth moments — (^qq) = 0.01 and (^^^j) = 0.07 — which shows breakdown of convergence 
of moment, at = 2fm~^ we can get (C.^^) = —0.37 and {^ao) ~ —0.13 which shows negative 
value but the convergence recovers. The moments and LCDAs of are plotted in Fig.[5]and 
right panel of Fig. El respectively. 

The instanton-involved twist-3 LCDAs calculated in this work present nontrivial properties, 
to some extent is unexpected. Although in principle non-positive-definite LCDAs are allowed, 
since here we lack of direct test of these LCDAs in exclusive processes we would like to give 
some tentative discussion on the possible ingredients which are not mentioned above and may 
have some impacts on our results. 

• multi-instantons 
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For simplicity we adopt single instanton approximation in our calculation, in fact there might 
be multi-instanton contribution to the correlation function. But the results in Ref. 



indicate 



that in singular gauge multi-instanton contribution to the pion correlator is coincident with that 
given by single-instanton approximation. While physical results should be gauge-independent 
thus we conjecture it might be reasonable to utilize single-instanton approximation in our 
calculation. On the other hand on can easily see it is very difficult to work out the instanton 
contribution in singular gauge for n 7^ 0. 

• instanton density 

It is obvious that the instanton density is crucial to obtain convergent results. Low instanton 
density is welcomed in our calculation. The density |fm~^ is widely used under single instanton 
approximation which give many reasonable results. In a way this indicates that low instanton 
density is consistent with single approximation. In other word it seems that high instanton 
density is questionable at single instanton approximation. Maybe it is the main reason that at 
high density the convergence is lost. 

• subleading Fock states 

From Eq.(2) one can see we use valence model to investigate twist-3 LCDAs. The nonlead- 
ing Fock states also may contribute twist-3 component via the mixing with other wave func- 
tions. [30]. This correction can be added by using the renormalization group, so it is less relative 
to this work. 



IV. CONCLUSIONS 

In the present work we have investigated the instanton effect by single instanton approx- 
imation on the twist-3 LCDAs of pion, /o(1370), J'rQ(1430) and ao(1450) from valence quark 
model within the framework of QCD moment sum rules. Results illustrate that the instanton- 
free twist-3 LCDAs are always positive-definite while the instanton-involved LCDAs show some 
nontrivial properties. We find that low instanton density is consistent with the method adopted 
in this work. Possible ingredients which might have impact on the results are briefiy discussed. 
Nonetheless we hope these LCDAs may be helpful to some heavy fiavored exclusive processes 
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since we conjecture the instanton density may play some role of a tuning parameter in deriving 
experiment favored results. 
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Appendix A: Vanishing of the pure zero-mode contribution to tensor moment sum 
rules 



We still work in four-dimension Euclidean space. Considering the following two-pint corre- 
lation function 

1 



d^xe''^'={0\TOn{x)O^{0)\0) 



2567r^ m*m2 
1 



ci^xe^^^ [ dpn{p)p^ I d^xa-^ 
J J + 



,2^3 



3/2 



with 



{x - Xq)^ + p2 

xtr [7a.7i^(1 - l5)(^af3lalp{l " Ts)] 



{x - Xo)^ + 



3/2 



(Al) 



On{x) = qio{x)aa^{iz ■ D /)"+^g2o(a:), O\0) = g2o(0)gio, 

where the integrations over instanton collective coordinates are explicitly. In fact it is enough 
to concentrate on the trace part in Eq. flAip 



tr 



T T T T 
/i u a p 



tr 



{5pu + ){8ap+iVbai. 



pi/^crp '^"Qapu^acrp 
= '^{Spu^ap ~ Spa^up + Spp^ua) ~ '^^pvap (^2) 

where the use of Eq. (123|) and Eq. (124|) have been made. Then combining above result with 
the trace over 7 matrix we have 



X 2tr 



''^^puap X 2tr 



(A3) 
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while in four- dimension Euclidean space 

1 

75 = ^£^,uapl^,7u7cT7p, 

Now one can see the whole trace part in Eq. (lAip vanishes thus there is no pure zero-mode 
contribution to tensor moment sum rules consequently the LCDA. 
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